r— l 

kO 

CO 


D D C 

l[jUZi2iZi]Di]E 


MAR  1979 

ETSEDTTlEi 
c 


Dynamics  of  Acoustic  Resonance  Scattering 
From  Spherical  Targets 
•• 

J.  George,  H.  Uberall,  A.  Farhan,  G. 

a 

Mezzorani,  A.  Nagl  and  K.  Sage 

Physics  Department,  The  Catholic  University 
of  America  - Washington,  D.  C.  20064 

and 

J.  D.  Murphy 

3709  Merlin  Way,  Annandale,  VA  22003 
Abstract 


31 

Q_» 

O 

C_3 


A resonance  theory  of  acoustic  scattering  from  penetrable 
bodies  has  recently  been  developed  (L.  Flax,  L.  R.  Dragonette, 
and  H.  Uberall,  J.  Acoust.  Soc.  Am.  63 , 1690  (1978))  which 
describes  the  features  of  the  scattering  amplitude  in  terms  of 
a background  amplitude  like  that  of  an  impenetrable  scatterer 
with  superimposed  resonances  at  the  known  eigenf requencies  of 
the  body.  in  the  present  study  (carried  out  for  the  example 
of  a spherical  gas  bubble  in  a fluid  which  in  spite  of  its 
simplicity  contains  all  the  essential  features  of  the  scattering 
process) , the  existence  of  these  resonances  is  explained  by  the 
phase-matching  of  surface  waves  on  the  scatter  whose  phase  and 
group  velocities  are  obtained  here.  In  addition,  we  study  the 
dynamics  of  the  interior  motions  inside  the  body  during  the 
scattering  process,  which  are  shown  to  undergo  resonances  in 
the  radial  and  angular  directions  except  for  the  radial  motions 
on  the  surface,  resulting  in  the  fact  that  the  modal  surface 
deformations  do  not  exhibit  any  particular  resonance  effects  at 
the  eigenf requencies  of  the  body 


a)  Now  also  at  the  Office  of  Marine  Technology,  Engineering 
Development  Laboratory,  NOAA,  Riverdale,  MD  20840 
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FOREWORD 

In  this  study,  we  investigate  the  dynamics  of  surface 
motion  of  a penetrable  spherical  body  subjected  to  an  inci- 
dent place  acoustic  wave.  The  results  should  be  valid  for 
objects  of  rather  general  internal  consistency  since  they 
are  based  on  quite  basic  physical  considerations,  but  they 
are  illustrated  here  by  the  example  of  a gas  bubble  in  a 
fluid.  It  is  shown  that  surface  waves  generated  in  the 
scattering  process  circumnavigate  the  target  sphere  and,  by 
a process  of  in-phase  "resonance  reinforcement",  excite  its 
multipole  resonances  at  those  eigenfrequencies  at  which  n + 1/2 
wavelenghts  (n  being  are  integer)  span  the  circumference. 

Phase  matching  is  achieved  here  due  to  the  fact  that  the  surface 
wave  suffers  a phase  jump  of  Tt  /2  at  each  of  the  two  "caustic 
points"  on  the  sphere  at  which  the  surface  waves  converge. 

Surface  motions  and  circumferential  waves  on  the  bubble  are  graphi- 
cally demonstrated,  both  on  and  off  resonance.  It  is  shown  that 
the  resonance  manifests  itself  by  a large-amplitude  angular  oscil- 
lation of  the  gas  contained  in  the  bubble. 


Introduction 


Recently,  some  progress  in  dynamical  studies  of  sound 
scattering  from  elastic  objects^  was  achieved  by  linking  the 
excitation  of  resonances  of  the  target  body  during  the  scattering 
process  to  a resonant  reinforcement  of  diffracted  "creeping 
waves"  in  their  multiple  circumnavigations  of  the  body.  In 
this  way,  the  occurence  of  acoustic  scattering  resonances,  and 
more  fundamentally,  the  very  existence  of  the  eigenvibrations 
of  an  elastic  body,  has  been  traced  back  to  a dynamic  origin 
rooted  in  the  phenomenon  of  circumferential  waves  that  exist 
on  the  surface  of  the  body,  and  partly  also  in  the  ambient 

medium.  Corresponding  resonances  appear  in  the  surface  waves 
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themselves  , and,  as  shown  by  Flax,  Dragonette  and  Uberall 

who  following  a suggestion  of  L.  Flax  used  the  methods  of 

nuclear  scattering  theory  for  this  purpose,  these  resonances 

become  manifest  also  in  the  individual  partial  waves  of  a 

normal  mode  expansion  of  the  field  where  they  interfere  with  a 

smooth  "potential  scattering"  background  characteristic  of  an 
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impenetrable  (rigid,  soft  ' , or  intermediate  ) target  object. 

Similar  results  were  also  observed  to  hold  in  the  scattering  of 
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elastic  waves  from  cavities 


The  link  between  creeping  waves  and  eigenvibrations  for  a 
spherical  target  takes  a somewhat  different,  and  more  interesting, 
form  than  for  the  cylindrical  bodies  discussed  earlier^" The 
eigenvibrations  and  related  creeping  waves  of  spherical  bodies  will 
thus  be  discussed  in  the  following,  considering  the  special  case 
of  gas-filled  bubbles  in  a fluid  medium.  Although  this  example 
constitutes  the  simplest  kind  of  scattering  problem,  it  exhibits 
all  the  basic  features  of  the  creeping-wave  resonances  to  be  dis- 
cussed here;  and  it  represents  an  interesting  problem  in  its  own 

right,  for  which  prior  to  our  recent  discussion'1''1",  only  scant  know- 
12 

ledge  existed  concerning  all  the  higher  multipole  resonances 

13 

(apart  from  the  we.  1-known  fundamental  "breathing-mode"  monopole 
resonance,  of  course) . 

It  will  be  shown  below  that  the  sound  scattering  phenomenon 
with  a spherical  target  bubble  can  be  resolved  into  the  creation 
of  a manifold  of  circumferential  waves  (represented  by  "Regge  poles" 
in  the  scattering  amplitude) , which  circumnavigate  the  target  with 
various  phase  velocities  that  depend  on  the  frequency.  Their  dis- 
persion curves  (as  well  as  their  trajectories  in  the  complex  mode 
number  plane)  may  be  obtained  from  the  calculated  progression  of 
the  modal  resonances  to  higher  frequencies  with  increasing  mode 
number,  as  studied  by  us  earlier^.  These  resonances  were  shown 


there  to  be  superimposed  on,  and  to  interfere  with,  a smooth 


background  amplitude  corresponding  to  the  scattering  from  a pressure- 
release  sphere,  for  the  case  of  an  air  bubble  in  water.  We  now  demon- 
strate  that  resonance  takes  place  in  the  nth  mode  at  any  such  frequency 
at  which  one  of  the  circumferential  waves  spans  the  circumference  of 
the  bubble  n + % times,  n being  an  integer.  It  is  shown  in  addition 
that  all  these  creeping  waves  converge  at  the  vertex  and  at  its  anti- 
pode on  the  sphere,  and  lose  a quarter-wavelength  (corresponding  to 
a phase  jump  of  fC/2)  upon  passage  through  each  one  of  these  two 
caustic  points.  One  full  circumnavigation  thus  leads  to  perfect 
phase  matching  when  the  above  condition  of  n + geometrical  wave- 
lengths equalling  the  circumference  is  satisfied,  and  hence  to  a 
resonant  reinforcement  of  the  creeping  wave  during  its  repeated 
circumnavigations.  This  phenomenon  forms  the  basis  for  the  resonant 
eigenvibrations  of  a gas  bubble  in  a fluid,  an  explanation  which  will 
also  hold  for  target  bodies  of  more  general  (e.g.  elastic,  layered, 
etc.)  consistency  and  for  the  excitation  of  their  eigenvibrations 
by  more  general,  non-acoustic  forces.  The  very  existence  of  the 
eigenvibrations  of  elastic  objects  may  then,  in  a sense,  be  referred 
back  to  the  existence  of  circumferential  waves  on  their  surface,  and 
to  the  resonant  reinforcements  of  the  latter  in  the  course  of  their 


repeated  circumnavigations. 


We  have  carried  out  numerical  calculations  in  order 


to  demonstrate  the  phase  matching  of  the  surface  waves. 

Further  numerical  calculations  were  done  for  study  of  the 

displacement  of  the  air  particles  in  the  bubble  on  and 

off  the  resonance.  It  is  found  that  at  resonance,  the  air 

3 

motion  in  the  bubble  is  enhanced  by  a factor  ***  10  , 

except  for  the  radial  component  of  motion  on  the  bubble 
surface.  This  results  in  the  fact  that  the  modal  bubble 
surface  deformations  do  not  exhibit  any  particular 
resonance  effects  at  the  eigenfrequencies  of  the  bubble. 
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I . Modal  Resonances  of  a Gas  Bubble  in  a Fluid 


The  scattering  of  sound  by  a fluid  (or  gaseous)  sphere  imbedded 
in  another  fluid  may  be  discussed  in  terms  of  a velocity  potential 
<f>  ( r)  , from  which  the  particle  velocity  field  v is  obtained  via 

V - - . (1) 

The  total  (incident  plus  scattered) field  is  given^1  by  the  normal- 
mode series 

o° 

d>  - fa  E (2tL+i)'tM*[^  (cos^)  (2) 

where  j is  due  to  the  incident  and  h ^ to  the  scattered  wave. 

Jn  n 

Here,  is  the  amplitude  of  the  incident  plane  wave,  and  k = 6l/c 

its  wave  number  is  the  ambient  fluid.  Satisfying  the  boundary  conditions 

(continuity  of  pressure  and  of  normal  particle  velocity)  at  the 

spherical  surface  r = a determines  the  scattering  coefficient  A of 

n 

the  nth  mode  as 


A = - 


ft  (X)  — Xy  (Xy)  >*.(■*') 


where  x = ka,  x^  = k^a.  with  k^  = w/c.^  being  the  wave  number  of  sound 
inside  the  bubble;  and  are  the  fluid  density  inside  and  outside 
the  bubble,  respectively. 
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Following  our  previous  study,  we  rewrite  Eq.  (3)  in  the  form 

l(sr-i) 


c*  (°)  g*  > h-l  A 

A « 2“  1 


AV 


Z Z'^-F 


(4) 


where 


'VC 

tx)  = 

(5a) 

X O0,  t=i,2 

(5b) 

7^  = 

/p*)  Xfe  )/^A  IXt)  • 

(5c) 

The  last  term  in  Eq.  (4)  corresponds  to  the  scattering  coefficient  for 

the  limit  of  a pressure-release  sphere,  where  p, /p  -»0  and  F -*co  . 

* b > w n 


Reference  11  further  showed  that  after  a Taylor  expansion,  Eq.  (4) 
becomes  approximately 


A ~ & 

'H. 


~ (i/Z) 


l"bi  X-\e->(i/z)rne 


+ -i-  Q.  Sin 


(6) 


where  x is  the  >^th  resonance  frequency  of  the  bubble  contained 
in  the  nth  normal  mode  of  the  scattering  amplitude  as  determined 
from  the  eigenvalue  equation 


Re  sO»  = Tt . 


(7) 


The  resonance  width  Pg  is  given  by 


^ — ( 2-/  P-ni  ) 2n  ^ , 


where 


6 = (d/dx) 

/ otC  x= 


Eq.  (6)  indicates  that  the  amplitude  consists  of  a smooth  background 
due  to  scattering  from  a pressure  release  sphere  (second  term) , super- 
imposed upon  which  there  is  a series  of  interfering  resonances  (first 
term) . This  is  demonstrated  in  Fig.  1 where  we  plot  the  quantities 
I U (*•)  I for  — TV  (-yv  S 3)  vs.  x 5 ka  up  to  x = 4 for 

an  air  bubble  in  water,  where 

(^)=  (i/k)(2x+l)  A^(co^) 

is  the  asymptotic  scattering  amplitude  (" form  function").  The  smooth 
background  is  clearly  visible,  and  the  interfering  resonance  peaks  are 
labeled  by  (n,£,).  The  tall  resonance  peak  (0,1)  represents  the  well- 
known^  giant  monopole  bubble  resonance.  The  resonance  frequencies, 
obtained  from  a numerical  solution  of  Eq.  (7)  , are  also  given  in 
Table  I up  to  x ^ 4. 
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II.  Circumferential  Waves  and  Their  Resonances 


14 

In  the  nuclear  scattering  theory  formalism  , one  usually  rewrites 
Eq.  (2)  as  <f?  = <f>^  + cpoc , where  the  amplitude  of  spherically  out- 
going waves  is 

CO 

- -f  6.-*“*  Y.  t 

'n=0 

with  the  "scattering  function" 


5L  - 2A„  + 1 . 

Due  to  Eq.  (6) , has  poles  in  the  lower  half  of  the  complex  x-plane, 
located  at 

x = ■ 

These  may,  however,  be  rewritten  in  the  form  of  poles  in  the  complex 
n-plane,  by  treating  n as  a continuous  variable  and  expanding 

+ }x*,*  +- 

( >0,  see  Tig.  i,  where  X^e  = (d/dtt)  ) , about  an  n-value  n. 

chosen  so  that  X 0 — X (the  given  incident  frequency).  The  S- 

r 

function,  which  reads 


d =r  C*  co)  ^ d'.ir  X-Ke  -U/2.)rMe 

* '*  2 fi)  qr  ^ Z_  ...  . i- 


I 


(note  that  the  various  terms  in  the  sum  are  accurate  only  in  the 
immediate  vicinity  of  their  resonance  peaks) , can  then  be  rewritten 


in  the  form 


^ (o)  ^ ***  (t/2.)  ^ (o)  ^ — 

ko  'n-*£-(i/i  )IZ/i  " 


i-o 


Ot-'VL, 


where 


f = P / x' 

' 'V  , 


14 


indicating  the  existence  of  "Regge  poles"  in  the  scattering  function 

A 

located  at  the  positions  yyi^  , where 


'yy.p  ~ 'yi.  + U/z)  fZ 


in  the  upper  half  of  the  complex  n-plane. 

These  Regge  poles,  and  their  residues,  are  conventionally  obtained 
by  performing  a Watson  transformation  on  Eq.  (10) , whereby  a series  of 


the  form 


oo 


14 


is  rewritten  as  a contour  integral: 


E t :p  rcosi*) —ti/i)  y l dv 

/Vi-o  ' J 1V-1  cos-rrv  * 

C being  a clockwise  contour  enclosing  the  positive  real  axis  in  the 
complex  v>  plane.  Applying  this  to  Eq.  (10)  with  an  appropriate 
expression  for  one  may  then  re-evaluate  the  contour  integral 


(15a) 


(15b) 


(15c) 


(16a) 
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at  the  Regge  poles  V = v ^ + T of  the  function Y , written  in  the 

V~T. 

form 


OO 


u - 1 


<p«> 

f 


J.o 

in  order  to  exhibit  the  denominator  in  Eq.  (15a)  . This  leads  to  the 
exact  expression 

7T  Kost5>) 


oo 


with 


fo.-t 


(2)  ^ 


C0S7T 


f'>V  i ) 


for  the  scattered  wave  of  Eq.  (10) , which  has  the  same  poles  and 
residues  in  the  variable  at  all  the  integers  n as  does  Eq.  (10) 

with  the  approximate  expression  of  Eq.  (15a) , and  we  see  that  our 
consideration  of  the  Watson  transformation  has  furnished  generalized 
resonance  denominators  which  interpolate  between  the  resonance  positions 
and  which  will  be  seen  to  provide  us  with  a physical  picture  of  the 
very  origin  of  the  resonances. 

14 

Such  a picture  may  be  obtained  by  using  the  asymptotic  form 
(valid  for  sin-$*.^  1)  , 

*/2_ 

ZP  (-cos  &)  = ( — cosfrv.i)  (Tr-^-^Trl 


fO) 


(16b) 


(17a) 


(17b) 


(18«q 
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in  Eq.  (17) ; even  for  small  values  of  n,  this  should  provide  an 
approximate  description  of  the  general  behavior  of  the  amplitude. 
One  may  then  split  the  cosine  in  Eq.  (18a)  into  two  exponentials, 
and  expand  the  denominator  of  Eq.  (17) , 

i/CosTT  2.  - 1 ! Sj'v,  TV^£  - 2 a.  **  ^2W+  1^‘*V 


'VHSO 


(18b) 


so  that 


OO 


• (19) 

^=0 


With  the  time  factor  exp(iw"t)  of  Eq.  (10),  this  shows  that  for  m = 0 

and  a given  value  of  , two  surface  waves  with  angular  propagation 
A 

constant  n + ^ propagate  in  the  counterclockwise  ( £ =-l)  and  in 
the  clockwise  ( £ = +1)  direction  around  a meridian  of  the  bubble, 
with  an  amplitude  spreading  factor  sin  that  describes  their  con- 

vergence at  the  north  and  south  poles  (i.e.  the  vertex  of  the  sphere 
as  seen  by  the  incident  wave,  and  its  antipode) . These  "tidal  waves" 
are  joined  by  other  such  waves  with  m > 0 that  already  have  encircled 

the  bubble  m times  previously  (of  course  with  larger  and  larger  atten- 

>\  . 

uation,  which  is  furnished  by  Im  n^  f 0) , corresponding  to  steady- 
state  scattering.  The  wavelength  of  these  surface  waves  is  obtained  as 


Af  (x)  ~ 2jCaX  j 1?e.  ('yLj+  k. ) 

■»  2.7T0C  / (*K y + i ) 


(20a) 


(20b) 
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if  Eq.  (15c)  is  used,  with  n determined  by  the  equation  X 

JL 

Due  to  radiation  damping,  these  waves  decay  as  exp(-'9yv>  ),  where 

0 ^ 
the  decay  angle  a7>  is  found  as 

f 


A (x>=  l/r 


■** 


-ix'  /r 

t ' ^ i 

if  Eq.  (15c)  is  used,  with  ^ and  X/  ^ given  by  Eqs.  (8a) 
and  after  Eq.  (13) , respectively.  Such  surface  waves  have  been  termed 
"creeping  waves"  by  Franz15,  although  their  discussion  goes  back  to 
van  der  Pol15  and  before. 

The  phase  velocity  of  the  sfcth  surface  wave  is  seen  from  Eq.  (19) 


(21a) 


(21b) 


to  be 


'q.(x)  = 


xxuo 


'Re  (*g  -*•  k. ) 


dispersion  curves  for  the  surface  waves  based  on  this  expression  will 
be  shown  below. 

Using  the  above  picture  of  circumnavigating  surface  waves  excited 

in  the  scattering  process,  one  arrives  at  a physical  explanation  for 

the  origin  of  the  resonances  themselves.  Equations  (20)  show  that  at 

A 

a physical  resonance  (Rert^  = integer  = n,  cf.,  e.g.,  Eq.  (15a)), 
the  wavelength  of  the  surface  wave  is  given  by 


-15- 


- 2T^/fw.  %.).  (23) 

A 

The  resonance  condition.  Re  n^  (x)  = n,  is  therefore  tantamount  to 
the  requirement  that  the  bubble  circumference  measure  an  odd  number  of 
half-wavelengths  of  the  surface  wave,  which  may  also  be  considered  the 
eigenvalue  equation  for  the  resonance  frequency  X^.  A resonance 
condition  of  this  kind  has  already  been  noted  by  Junger  and  Feit^7. 

Equation  (23)  actually  leads  to  phase  matching  of  the  surface  wave  after 
each  additional  curcumnavigation,  in  spite  of  the  fact  that  a half- 
integer number  of  wavelengths  fit  around  the  bubble  circumference. 

This  becomes  clear  from  inspection  of  Eq.  (19) , which  reveals  that 
owing  to  the  factors  exp  (-i  £1T /4)  , a phase  jump  corresponding  to 
exp(-  lit  /2)  takes  place  every  time  the  surface  wave  passes  through  the 
convergence  zones  at  the  north  and  south  pole  of  the  spherical  bubble 

($  = 0/  7T  ) . During  such  a passage,  the  continuing  description  of 

18 

a surface  wave  encircling  the  bubble  in  a given  sense  passes  from 
one  exponential  in  Eq.  (19)  to  the  other  one,  the  relative  difference 
in  the  factors  exp(-  elf  /4)  providing  the  mentioned  phase  jump. 

Since  the  latter  amounts  to  a backVAfd-  shift  of  the  wave  by  an 

extra  % + \ wavelengths  (at  each  pole)  = % wavelength  during  one  complete 

circumnavigation,  the  condition  2 Tf  a.  = ( re  + *5)  A.  actually  results  in 
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an  exact  phase  matching  of  the  surface  wave  with  itself  after  each 
complete  circuit,  which  explains  physically  the  appearance  of  a resonance 
as  a resonant  reinforcement  of  the  tidal  waves  in  the  course  of  their 
repeated  circumnavigations.  If  unattenuated,  such  a reinforcement 
would  cause  infinite  amplitudes.  The  radiative  attenuation  of  the 
surface  wave  leads  to  finite  amplitudes  and  a finite  resonance  width. 

As  mentioned,  the  north  and  south  pole  represent  caustic  points 
for  the  surface  waves  on  the  bubble.  A corresponding  phase  jump  of 
-7 X /4  for  a wave  passing  through  a caustic  is  well-known  for  other 
examples  in  acoustics,  e.g.  sound  propagation  through  the  inhomogeneous 
oceanic  medium"^. 

The  phase  matching  condition  may  be  viewed  somewhat  differently, 

namely  as  a coincidence  condition  for  phase  velocities.  The  vibrational 

modes  continued  in  Eq.  (10) , being  standing  circumferential  waves 

around  the  bubble  circumference,  may  each  be  decomposed  into  a pair 

of  modal  waves  propagating  in  opposite  directions  with  phase  velocities 

20 

€.  , using  the  separation 


~ Q,  k |Pv(-Cojv^+  Tp Qv 

= [z/tC  (V  + i)  \ ^ tV*  ^ J 

as  shown  in  Fig.  2 for  the  quadrupole  (n  = 2) . The  acoustic  scattering 
process  generates  a series  of  surface  waves  with  speeds  c . (x)  around 
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the  bubble  given  by  Eq.  (22),  while  using  Eqs.  (10)  and  (18),  the  modal 
speeds  are  given  by 


/C  = XX.  CO  /C'W  Z.  ) • 

/iV  # 


(25) 


From  Eq.  (22) , it  is  seen  that  a resonance  occurs  at  that  frequency 
at  which  the  speed  of  the  th  surface  wave  launched  during  the 
scattering  process  coincides  with  the  phase  velocity  of  the  modal  wave 
corresponding  to  the  nth  natural  multipole  vibrational  mode  of  the 
bubble.  Again,  this  coincidence  condition  for  phase  velocities  may 
be  viewed  as  an  eigenvalue  condition  for 


III . Numerical  Results 

In  the  following,  we  shall  obtain  some  numerical  results  of  the 
properties  of  surface  waves  as  they  are  excited  by  plane  acoustic 
waves  incident  on  an  air  bubble  in  water. 

1. ) Phase  and  group  velocities. 

The  phase  velocity  of  the  /£th  surface  wave  is  given  by  Eq.  (22) . 

In  Fig.  3a,  we  plot  the  ratios  c.  /c  as  a function  of  x = ka  for  ^=1 

through  14.  The  ensuing  dispersion  curves,  which  are  not  unlike  those 

21 

for  the  Lamb  waves  or  a fluid  layer  in  a vacuum  , come  in  from  infinity 
at  some  low-frequency  cutoff  (except  for  the  X = 0 wave),  and  at  high 
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frequencies  they  allappear  to  approach  the  common  value 

c^/c  = 0.22,  where  c^  is  the  sound  velocity  in  the  fluid 

contained  in  the  bubble  (here,  air) . 

If  we  introduce  a trace  velocity  c = c/sin  Ot  of 

the  incident  plane  wave  on  the  surface  of  the  sphere, 

where  is  the  angle  of  the  surface  normal  with  the 

direction  of  incidence,  then  Cremer 1 s "coincidence 

between  the  incident-wave  trace 

velocity  and  the  surface  wave  speed  determines  the  angle 

CX  at  which  the  surface  wave  is  launched.  The  infinite 

values  of  c^  at  the  cutoffs  in  Fig.  3a.  thus  correspond 

to  a launching  angle  & = 0.  At  those  frequencies  where 
/C ^ < /C  • the  launching  presumably  then  occurs  at 
lancing  incidence,  c<  = 9o°  , if  an  excitation  of 

these  waves  does  take  place  at  all. 

The  group  velocity  of  the  jtih  surface  wave  is  given 


22 

condition"  c^  = c 
tr  Z 


by 


XX. 


d * 


oo  = — 

**  xi  Tse  (y^  * 

In  Fig  , we  plot  the  ratios  as  a function  of 

x £ for  ^ =1  through  14.  They  were  obtained  from 

our  calculated  resonance  frequencies  which  only  covered  a 
limited  region,  so  that  some  portions  of  the  graph  remained 
open. 

2 . ) Circumferential  Waves 

We  shall  here  illustrate  the  propagation  of  the  previously 


(26) 


I 


m 
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discussed  circumferential  waves  around  the  bubble  during  the 
scattering  process,  exhibiting  the  -KjZ  phase  jump  at  the 
poles  of  the  sphere  and  the  phase  matching  at  resonance. 

Note  that  the  field  of  Eq.  (19)  actually  describes  the 

total  field  <j>  ( = vrL’hCj>  outside  the  bubble,  since 

(P does  not  contain  the  S function  with  its  Regge 
poles,  hence  furnishes  no  contribution  in  the  residue 
evaluation  of  the  Watson  integral. 

After  expanding  the  denominator  of  Eq.  (17a)  as  in  Eq. 

(18i>)  , we  obtain  using  Eq.  (24): 

©O  <*3 

with 


(27a) 


£,tn- 0 


L (p,^)  -=  It  l i)7rL^p^ (-**&)  + (271 


or,  asymptotically: 


\1'tpW  (2».l)«rt,  ± J (27c) 


,£  Z 


where  0(  = Tl  - , and  P=2m7T'  + ot  is  the  continuosly 

increasing  phase  angle  around  the  sphere,  cf.  footnote  18. 
For  a given  label  £ characterizing  the  order  of  the 
surface  wave,  this  describes  two  waves  encircling  the  bubbfe 
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in  opposite  directions  with  angular  propagation  constant 

-L  . We  shall  represent  graphically  one  of  these 


A 

'VL/i  + 


waves  by  the  surface  deformation  which  it  causes.  The 
deformation  is  determined  by  the  deformed  radius^ 

= -a  + "Re  (*■,$), 

and  the  radial  displacement  is  given  exactly  by 


Re 


■ l ($r  F.  (y^l)  I 


T/nzO 


or  asymptotically  ( for  /visin'!/*  ^ 1 ) by 


rVi 


Re  u«f*  («.*)  = -I.  F0  /^  (,)|  (*.«' 

* Cos  ^ ± (0*.+  4 )L°*  ± + ^ ^ \ > 


with  an  amplitude  factor 


= [ Ux)1'7<l]  ji.  , 


where  due  to  the  narrowness  of  the  resonances  (Fig.  1)  and 
the  ensuing  weak  attenuation,  we  were  able  to  approximate 
= 71 q for  the  case  of  the  air  bubble  in  water. 

In  Fig.  4 and  5,  we  plot  the  amplitude  of  the  surface 
deformation  caused  by  circumferential  waves  as  given  by  Eq. 


(28a) 


(28b) 


(28c) 


( 28d) 


(28^  )»  tor  various  cases  On  and  off  resonance.  The  plots 
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are  done  on  an  arbitrary  scale,  only  the  wave  which 
progresses  in  the  clockwise  sense  is  shown,  and  ( x ) \ 

is  approximated  by  a constant  in  the  vicinity  of  each 
resonance.  In  addition,  we  multiply  the  displacement 

U)  , ft  h 

amplitude  Reu^  by  ( st-nir*  ) in  order  to  remove  the 
corresponding  singularity  which  is  merely  caused  by  the 
geometrical  convergence  of  the  waves  at  the  poles  of  the 
sphere. 

Fig.  4a  shows  a surface  wave  with  n = 4.7,  so  that  we 
are  closfe  to  the  resonance  case  n^  = 5 at  which  5.5  wave- 
lengths span  a meridian  of  the  sphere.  It  is  evident  that 


at  each  of  the  poles  at 


$ =o°  and  180°  a phase  jump 


by  a quarter  wavelength  takes  place,  and  also  that  after  one 
circumnavigation  the  wave  does  not  close  into  itself,  so  that 
no  resonant  buildup  of  the  amplitude  will  result.  The  solid 
curve  has  been  obtained  using  the  asymptotic  form  Eq.  (28c  ) 
which  is  quite  adequate  for  this  value  of  /rt  , and  the 
dashed  curve  (shown  here  only  in  the  vicinity  of  & = ISO0  ) 
was  obtained  from  the  exact  Eq.  (28  d)  in  which  and 

Q with  non-integer  index  were  calculated  numerically 

^ 22 
in  terms  of  hypergeometric  series.  It  is  seen  that  in  tie 

exact  calculation,  the  phase  jump  is  not  abrupt  but  takes 

place  in  a smooth  fashion. 


P 
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In  Fig.  4b,  the  same  quantities  are  plotted  for  n = 5.0, 

X 

so  that  the  surface  wave  closes  into  itself  and  a resonance 
results.  We  see  that  at  & = 180*  , the  exact  displacement 
vanishes;  this  takes  place  in  such  a fashion  that  if  the  geometri- 
cal amplitude  factor  (sin  oh  ) 2 is  retained,  the  resonating 

amplitude  remains  finite  at  & = 180°  . This  is  shown  by  the 
light  dashed  line  on  the  extreme  right  of  the  figure,  while  the 
light  solid  line  represents  the  asymptotic  calculation,  in  which 
the  factor  ( sfn'&  ) 2 causes  the  geometrical  divergence  of  the 

amplitude  which  is  evident  here. 

Analogous  results  are  shown  in  Fig.  5 near  and  at  the 
= 1 resonance  where  1.5  wavelengths  fit  over  the  meridian 
of  the  sphere.  Fig.  5a  shows  the  surface  wave  below  the 
resonance,  at  n^  = 0.7,  and  Fig.  5 \>  sit  the  resonance, 

= 1.0.  Again  the  solid  line  results  from  the  asymptotic  and 
the  dashed  line  from  the  exact  calculation. 

3 , ) Displacements  of  bubble  surface  and  interior  fluid 

The  particle  displacement  in  the  sound  field  is  related  to 
the  velocity  field  by  AJL  = , so  that 

XL  * (<-/co)3)-  = (iu>)  ^ , (29) 

see  Eq.  (1) . The  change  of  shape  of  the  bubble  surface  is 
governed  by  the  radial  displacement  u^,  which  is  continuous 


V 


r = ' 
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1 
1 

across  the  surface.  (The  angular  displacement  , -to  first 

order,  only  moves  the  air  particles  along  the  surface  without 

shape  change.  It  is  discontinues  across  the  surface,  since 

the  air  can  slip  relative  to  the  water  particles) . j 

Using  the  scalar  potential  within  the  bubble,'1'1  i 

co 

<&,  = 4>a  e"i“t  Y.  (V)A  («s$>,  (30a) 

'n  = o 

where 

we  find  for  the  displaced  radius  R of  the  bubble  in  its  extreme 

'rv 

positions,  if  only  the  nth  mode  is  excited: 

'R^  (-a,^)  = -a-  - (<£0/oc6o)  (2^  + i.)  (H)'?* (<*&).  (31) 

This  has  the  typical  shape  of  a spherical  multipole  of  nth 
order  ; it  is  plotted  (choosing  the  positive  sign)  in  Fig.  6 for 
n = 0 (monopole  or  breathing  mode,  a radially  symmetric 

in-and-out  pulsating  motion),  n = 1 (dipole  or,  to  first  order, 
rigid  back-and- forth  motion) , n=2  (quadrupole,  or  alternation 
between  prolate  and  oblate  ellipsoidal  shape),  and  n-=  3 
(octupole,  or  pear-shaped  deformation).  The  plots  in  Fig.  6k 
are  made  at  the  fundamental  resonance  frequency  of  each  mode11 
(shown  in  Table  I)  , and  in  addition  at  3 % below  and  above 

each  resonance  frequency,  respectivley . The  scale  is  arbitrary 


but  is  the  same  for  a given  value  of  n. 


It  is  seen  that  no 
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pronounced  (i.e.  factors  '“fw/f'i  ~I0J  ) resonant  enhancement 
takes  place  in  the  radial  displacements  of  the  surface,  while 
some  moderate  (factors  ~ 10)  enhancement  does  seem  to  occur. 

Mathematically,  this  arises  as  follows.  As  shown  in  Ref. 

11  (Eq.  35) , one  has  at  resonance: 

jj.  (xfc)  = Or  , (32) 

so  that  in  Eq.  (31)  will  be  of  order  , while  off 

resonacnce,  it  is  of  order  unity.  However,  only  the  combination 
(Xj,)  appears  in  Eq.  (31)  so  that  ■$  ) remains  of 

order  unity  at  resonance  and  does  not  show  any  significant 
resonance  effects.  The  physical  reason  for  this  lies  in  the 
boundary  condition  of  continuity  of  the  radial  displacement, 
linking  the  interior  radial  motion  to  that  of  the  external 
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fluid  on  the  boundary.  Since  only  the  irferior  fluid  reson  tes 
at  the  resonance  frequencies  of  the  bubble,  this  linkage 
prohibits  a pronunced  resonance  motion  of  the  shape  of  the 
bubble  surface. 

VJe  now  proceed  to  plotting  the  motion  of  the  air  particles 
in  the  interior  of  the  bubble,  and  shall  find  that  it  is  here 
where  the  pronounced  ( order  10  ) resonance  effects 

occur.  This  shows  clearly  that  the  bubble  resonances  are 
resonances  of  the  interior  fluid.  First,  we  consider  the  radial 
particle  displacements  inside  the  bubble  which  may  be  described 
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by  the  deformation  of  an  originally  spherical  surface  of 
radius  r;  this  will  be  changed  in  the  extreme  position  according 

40 -r  + (2*i+J.)  IB^I  (33) 

In  Fig.  7,  we  plot  the  quantity  R (r,o)  - T* , i.e.  the 
displacement  at  iJ>=  0 (where  it  assumes  its  maximun  value  since 
*1?  (1)  =1)  as  a function  of  r.  This  done  for  the  cases  ) 

equal  to  (o,l)  and  (0,2),  i.e.  monopole  fundamental  and  first 
overtone  (Fig.  7a),  further  (1,1)  and^l,2),  i.e.  dipole 
fundamental  and  first  overtone  (Fig.  llo  ),  and  finally  (2,1) 
and  (2,2),  i.e.  quadrupole  fundamental  and  first  overtone 
(Fig.  7c) . As  before,  the  displacement  is  drawn  at  the 
corresponding  resonance  frequency  of  Table  I,  and  at  a 
frequency  3%  below  (-)  and  above  (+)  resonance. 

The  following  information  is  obtained  from  this  figure: 

(a)  A resonance  enhancement  takes  place  in  the  radial 

in1 

displacement  inside  the  bubble,  which  is  of  order  ^ w /P  b ^ 

2 

(or  at  least  10  ) while  no  comparable  enhancement  occurs  at  the 
surface.  The  fundamental  monopole  (0,1)  is  an  exception  to 
this. 

(b)  Not  counting  any  nodes  in  the  low -amplitude  portion 

near  the  surface,  the  number  of  radial  nodes  N in  the  radial 

r 

displacement  is  given  by  N^_  = Ji  - 1 (the  monopole  n = 0 being 
an  exception) . Generally,  Nr  is  determined  by  the  behavior  of 
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as  a function  or  r,  for  at  the  appropriate 

resonance  frequency. 

(c)  For  higher  multipolarity  n,  the  pattern  of  large  radial 
displacements  shifts  closer  to  the  surface. 

(d)  Due  to  the  narrowness  of  the  (n,  £ ) 4=  ( 0,1) 

resonances,  a - 3%  frequency  shift  brings  us  well  away  from  the 

resonance  while  for  the  monopole  fundamental  ( 0 ,1),  it  still 

leaves  us  high  up  on  its  shoulder, 
ed 

The  mention/resonant  radial  enhancement  occurs  mathematically 
because  the  factor  in  Eq.  (33),  which  shows  the  mentioned 

enhancement  as  the  frequency  passes  through  the  resonance 
frequency,  is  now  accompanied  by  ^(^rj  rather  than  by 
^ ^ and  only  the  latter  factor  has  a compensating  minimum 

at  resonance  as  mentioned  above. 

In  Fig.  8,  we  present  some  graphics  of  the  full  angular 
pattern  of  radial  displacements  directly  at  resonance,  plotted 
at  originally  equidistant  radial  intervals.  This  is  shown  for 
the  case  (n ,Ji)  = ( 0^1)  and  (0,2)  in  Fig.  8a  (monopole  funda- 
mental and  its  first  overtone),  for  (1;1)  and  (1;2)  in  Fig.  8t 
(dipole  fundamental  and  its  first  overtone),  and  for  (2,1)  and 
(2,2)  in  Fig.  8/C  (quadrupole  fundamental  and  its  first  overtone), 
indicating  the  uniform  radial  shifts  that  take  plance  in  the 
fundamental  mode  and  the  non-uniform  shifts  in  the  overtones. 
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We  finally  consider  angular  displacements  of  the  bubble 
contents,  and  the  resonance  effects  that  occur  in  fhem.  Prom 
Egs.  (29)  and  (30),  the  extreme  values  of  the  modal  angular 
displacements  of  the  air  particles  inside  the  bubble  are 
obtained  as 

(Y]  & [^0/r-uo)  (2-u+O  IB^l  ^ ('sr)  SlVS'P,, 

The  argument  after  Eq.  (32)  now  indicates  that  will 

also  show  a resonance  enhancement  as  the  frequency  passes 
through  the  resonance  frequency  since  the  coefficient  B ^ 
does;  the  latter  is  here  accompanied  in  Eq.  (34)  by  the 
factor  j ^ ( k^r  ) which  is  of  order  unity  throughout,  so 

that  unlike  R ^ ) , the  angular  particle  motion  shows  a 

resonance  maximum  at  the  resonance  frequency  even  at  the 
bubble  surface. 

In  Fig.  9,  we  plot  (y-($  ) indicated  by  arrows, 

for  a number  of  appropriate  values  of  ^ and  at  equidistant 
radial  intervals,  for  the  modes  n = 1 and  £ (no  angular 
motion  takes  place  for  the  monopole  n = 0 which  constitutes  a 
radial  in-and-out  pulsation)  st  the  fundamental  resonance 


(34) 


frequency  of  each  mode 


and  at  the  first  overtone 


28 


s 


frequency.  All  scales  are  relative  here. 

It  is  seen  that  the  dipole  has  angular  modes  at  0 ° 

and  180°  , and  maxima,  at  and  270°  ; the  quadrupole 

has  nodes  at  = 0°  , 90°  , 180°  and  270°  and  maxim  at 

'S'  = 45°  , 135°  , 225°  and  315°.  The  octupole  has  nodes 

at  ■9'=  0°  , 60 ^ etc.  and  maxima,  at  30°  , 90°  etc. 

In  Fig.  10,  we  plot  the  angular  displacement  of  the  modes 
n=l,2  and  3 at  the  angle  of  its  maximum  ys.  the  radial 
coordinate  ; in  Fig.  10  a for  the  fundamentals  of  the  n = 1,2 
and  3 modes,  in  Fig.  101a  for  the  first  overtones.  We  also 
mention  that  as  in  the  case  of  the  radial  displacements,  the 
values  of  4JL  (f,  $ ) drop  by  values  of  the  order  10  2 when 

we  move  the  frequency  away  by  i 3%  from  the  resonance 
frequency;  this  effect  is  not  graphically  represented  here. 

Figs.  9 and  10  together  with  what  was  just  said,  indicate 
that  as  stated  above,  the  angular  displacements  exhibit 
resonance  effects  throughout  the  interior  of  the  bubble,  right 
up  its  surface.  At  that  place,  the  air  motion  slips  against 
the  motion  of  the  water  particles  which  does  not  show  any 
resonance  effects. 


Our  foregoing  discussion  of  the  radial  and  angular  displace- 


ments  of  the  air  content  of  a bubble  indicates  that  a resonance 
manifests  itself  here  by  a violent  motion  of  the  air  particles 
inside  the  bubble,  rather  than  by  a resonant  deformation  of  the 
bubble  surface.  It  can  be  expected  that  similar  results  hold 
for  the  interior  motion  of  the  elastic  medium  in.  the  case  of 
resonant  sound  scattering  from  a solid  object. 
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IV.  Conclusions 

We  have  analyzed  the  dynamics  of  the  resonant  behavior  of 
a penetrable  spherical  body  which  is  subjected  to  an  incident 
plane  acoustic  wave,  in  terms  of  circumferential  waves  as 
described  by  Regge  pole  theory.  Considering  the  air-filled 
spherical  bubble  immersed  in  water  as  an  example,  it  was 
shown  that  an  infinite  set  of  Regge  poles  labeled  by 
2 = 0,1,2...,  moving  along  trajectories  in  the  complex 
n-plane  with  increasing  frequency,  generate  in  succession  the 
fundamental  (£  = 0)  and  overtone  ( J.  ^ 1)  modal  resonances  of 
succesively  higher  multipolarity  n every  time  a pole  moves 
past  the  integer  value  n. 

We  have  shown  that  the  Regge  poles  cause  the  appareance 
of  circumferential  surface  (or  tidal)  waves  on  the  sphere 
that  are  generated  due  to  the  action  of  the  impinging  sound, 
but  which  could  equally  well  have  been  created  by  a more 
general  external  force  since  they  are  a property  of  the  scattering 
object  and  not  of  the  particular  excitation  mechanism.  The 
occurence  of  the  resonances  themselves  was  then  recongnized 


to  take  place  at  frequency  at  which  an  integer-pLus-half  number 
of  wavelengths  (the  irteger  being  equal  to  the  multipolarity  n 
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of  the  collective  bubble  vibration)  could  fit  over  the 
sphere's  circumference,  so  that  the  phase  matching  of  the 
repeatedly  circumnavigating  surface  waves  (taking  into 
account  a phase  jump  of  — 7i/2  of  the  surface  waves  when 
passing  through  each  of  the  convergence  zones  at  the  north 
and  south  pole)  will  lead  to  their  resonant  reinforcement. 

Graphical  illustrations  were  presented  for  the 
propagation  of  these  surface  waves,  as  well  as  of  the  modal 
surface  deformation  and  of  the  gas  motion  inside  the  bubble. 
The  latter  manifestly  exhibit  the  resonance  effects,  while 
the  surface  deformations  do  not,  on  account  of  the  boundary 
conditions . 

The  here-analyzed  dynamic  effects  around  spherical 

gas  bubbles  in  water  are  expected  to  be  of  very  general 

validity,  and  to  be  applicable  to  spherical  scatterers  of 

much  more  general  nature  (e.g.,  elastic  bodies  in  fluids, 

23 

fluid-filled  cavities  or  solid  inclusions  in  an  elastic 

24  25 

medium  etc.)  In  fact,  we  have  been  able  to  describe  in 

very  similar  terms  the  phenomenon  of  the  giant  collective 

resonances  of  atomic  nuclei  which,  on  the  basis  of  a 


hydrodynamic a 1 model  of  the  nucleus,  may  also  be  thought 
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of  being  caused  by  the  phase  matching  of  circumferential 
"nuclear  tidal  waves"  generated  by  various  kinds  of  nuclear 
reactions . 
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Table  I 

Resonance  frequencies  of  the  fundamental  U = 1) 

and  overtones  (.1  = 2,  3,  ...)  of  the  lowest  multipole  resonances 
(n  = 0 ...  ^1)  of  an  air-bubble  in  water. 


multipole 

Hi 

1 

2 

3 

4 

5 

6 

monopole 

0 

0.01391 

1.0161 

1.7467 

2.4658 

3.1792 

3.8938 

dipole 

1 

0.4712 

1.3433 

2.0817 

2.80504 

3.5229 

— 

quadrupole 

2 

0.7555 

1.6485 

2.4000 

3.1311 

3.8539 

— 

octupole 

3 

1.0214 

1.9411 

2.7074 

3.4478 

— 

— 
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Figure  Cap  Hons 

Fig.  1.  Moduli  of  individual  partial-wave  backscattering 

amplitudes  |f  (>.Tr)|  (-a  = 0 3)  plotted  VS. 

j for  an  air  bubble  in  water,  showing  resonances 
superimposed  upon  a smooth  background. 

Fig.  2.  The  nth  (here  sketched  as  2nd)  normal  mode  vibration 
decomposed  in  a pair  of  traveling  waves  with  phase 
velocities  C and  surface  waves  labeled  by 
^ , with  phase  velocities  c^  • 

Fig.  3.  (a)  Phase  velocities  of  surface  waves  Jl  - 1 

through  14  as  a function  of  /cd.  (dispersion  curves) 

relative  to  sound  velocity  in  water,  /C  , for  an 
air-filled  bubble  in  water.  (b)  Group  velocities 
C/  / /C  of  the  same  surface  waves. 

Fig.  4.  (a)  Surface  wave  amplitude  on  a sphere  slightly  below 

a resonance  ( /Yl ^ = 4.7),  with  no  phase  matching. 

Quarter-wavelength  phase  jumps  a-t  the  poles  are  evident, 
(b)  Surface  wave  at  the  = 5 resonance,  with  phase 
matching.  Solid  line:  asymptotic  calculation?  dashed 
line:  exact  calculation.  The  heavy  curves  disregard, 

and  the  light  curves  reta-in,  the  convergence  factor 


( Si'n  & ) 


-1/2. 


Fig.  5.  Same  as  Fig.  4,  for  ( a ) 'VL.  = 0.7  and  (b)  'Yip  = 1, 
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Fig.  6.  Surface  displacements  of  an  air  bubble  in  water,  set  into 
oscillation  by  impinging  sound,  for  the  modes  (a)  n = 0 
(monopole  or  breathing  mode)  , (b)  n = 1 (dipole)  , (c) 

n = 2 (quadrupole)  and  (d)  n = 3 (octupole) . Plots 
are  at:  3%  below  the  fundamental  resonance  frequency 

given  in  Table  I;  at  the  resonance  frequency,  and 

at  3%  above  the  resonance  frequency  (fundamentals 
only) . Scales  are  arbitrary,  but  are  the  same  for  a 
given  value  of  n. 

Fig.  7.  Radial  displacements  inside  the  bubble  as  functions 

of  r,  at  the  appropriate  resonance  frequency  as  well  as 
3%  below  ( - ) and  above  ( + ) the  latter:  (a)  for 
monopole  fundamental  (0,1)  and  its  first  overtone  (0,2), 

( 'b  ) for  dipole  fundamental  (1,1)  and  its  first 
overtone  (1,2),  and  ( c ) for  quadrupole  fundamental 
(2,1)  and  its  first  overtone  (2  ).  Scales  are  arbitrary, 
but  are  the  same  for  a given  (n , ji)  . 

Fig.  8.  Resonant  displacement  of  originally  equidistant  spherical 
surfaces  inside  the  bubble,  plotted  at  the  resonance 
frequency  for  (a)  the  monopole  fundamental  (n, C)  = (0;1) 
and  first  overtone  (0;2),  (b)  the  dipole  fundamental 

('Tt/  2 ) = (1/1)  and  first  overtone  (1,2),  and  (c)  the 
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quadrupole  fundamental  (2;1)  and  first  overtone  (2,2)  - 
arbitrary  scales. 

Fig.  9.  Angular  displacements  u^(r,/\)')  indicated  by  arrows  (ar- 
bitrary scale)  for  n = 1 and  2 at  the  fundamental  re- 


sonance frequency,  i.e.(l,l)  and  (2,1),  and  at  the  first 
overtone,  i.  e.  (2,1)  and  (2,2). 

Fig.  10.  Angular  displacements  u^,(r,^)  of  the  modes  n = 1,  2 and 
3 vs.  r/a,  each  taken  at  the  angle  ^ at  which  maximum 
displacement  occurs. 


